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SUMMARY 

In  this  continuing  program  a  small  disposable  laser  altimeter  developed 
earlier  in  the  program  was  used  to  measure  the  reflected  signal  from  the 
sea  surface  as  represented  by  the  waters  beneath  the  Golden  Gate  Bridge. 
The  work  revealed  a  number  of  new  features  of  the  reflection  not  previously 
reported  elsewhere.  The  bridge  was  used  as  a  site  both  because  the  cost  of 
helicopter,  airplane  or  blimp  platforms  would  have  been  prohibitive  and 
because  sampling  sequences  at  precise  angles  from  the  vertical  would  have 
been  difficult  to  achieve  from  such  platforms.  Thirty-second  samples  of 
the  return  signal  were  taken  at  a  one  kilohertz  rate  at  each  of  a  number  of 
angles  from  the  vertical  in  the  upwind  and  crosswind  directions,  during 
twelve  separate  occasions  on  the  bridge.  The  probability  distribution  for 
the  frequency  of  occurrence  of  return  signals  as  a  function  of  signal  amp- 
litude snowed  interesting  new  structure.  In  the  presence  of  swell  superim- 
posed on  a  shorter  wavelength  pattern,  periods  of  zero  return  signal  occur, 
even  for  vertical  incidence.  This  is  apparently  due  to  tipping  of  the 
reflection  scatter-cone  of  the  smaller  waves  far  enough  from  the  vertical 
that  no  light  is  reflected  vertically.  The  probability  distribution  of  the 
nonzero  signals  in  this  situation  has  a  shape  that  can  be  fitted  with  an 
antilog-normal  distribution.  This  is  skewed  toward  higher  signal  than  that 
provided  by  a  normal  (Gaussian)  distribution.  When  the  incidence  angle  is 
displaced  from  the  vertical,  the  distribution  shape  stays  the  same  but  the 
zero  reflections  become  more  frequent.  The  average  signal,  including  the 
zeroes,  falls  off  with  angle.  The  variation  of  the  average  signal  as  a 
function  of  angle  is  well  fitted  with  a  Gram-Charlier  distribution,  the 
same  as  found  earlier  by  other  observers  using  photographic  techniques 
which  did  not  reveal  these  details  of  the  structure.  The  Gram-Charlier 
distribution  is  somewhat  more  peaked  than  Gaussian  in  both  the  crosswind 
and  up/downwind  direction.  In  the  up/downwind  direction  it  is  slightly 
skewed  upwind.  For  the  simpler  wave  pattern  produced  by  a  single  long 
sustained  wind  direction,  the  probability  distribution  of  the  signal  ampli- 
tude is  log-normal  and  no  periods  of  zero  reflection  occur.  For  this  situ- 
ation the  distribution  shifts  from  lognormal  toward  .exponential  as  the 
angle  from  the  vertical  is  increased.  At  large  angle  the  distribution  is 
well  fitted  by  an  exponential  curve.  For  intermediate  surface  states,  with 
a  mixture  of  waves  of  different  amplitudes  and  different  directions,  but  in 
the  absence  of  large  swell,  the  distribution  is  usually  intermediate  be- 
tween lognormal  and  antilog-normal,  being  normal  in  some  cases.  The  over- 
all behavior  of  the  reflection  is  more  favorable  for  disposable  laser 
altimeters  than  previously  believed  because  the  skew  toward  higher  ampli- 
tudes gives  frequent  large  signals,  higher  amplitudes.  Also  for  an  altime- 
ter which  may  be  swinging  from  a  parachute  or  balloon,  the  return  at  angles 
from  the  vertical  remains  high.  The  presence  of  occasional  zero  return 
signal  does  degrade  the  accuracy  of  altitude  somewhat  for  a  descending 
altimeter,  but  the  enhanced  signal  due  to  the  distribution  skew  improves 
performance  at  larger  altitudes. 
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I.  INTRODUCTION 

The  work  reported  here  provides  detailed  information  on  the  nature  of  the 
reflection  from  ocean-like  waves  for  narrow  beam  laser  altimeters  or  simi- 
lar devices.  The  work  was  a  continuation  of  a  project1  in  which  a  small, 
light,  rugged,  and  inexpensive  disposable  laser  altimeter  was  developed  for 
parachute  or  balloon-borne  use  over  the  ocean.  In  such  a  device,  laser 
power,  detector  sensitivity  and  signal-to-noise  considerations  dictate  use 
of  as  narrow  a  laser  beam  and  detector  field  of  view  as  possible.  In  the 
course  of  field  tests  it  became  apparent  that  narrow  laser  beam  reflection 
from  ocean-like  water  waves  was  different  from  that  previously  reported  for 
wide  angle  laser  or  radar  beam  reflection.  The  reflection  was  found  here 
to  be  more  favorable  for  narrow  beam  laser  altimeters  than  for  the  log- 
normal  probability  distribution  of  intensity  that  is  generally  assumed. 
The  results  indicate  that  a  log-normal  distribution  is  applicable  to  wide 
angle  beams,  or  for  narrow  beams  with  short  wavelength  water  waves  in  the 
absence  of  long-wavelength  swell.  For  narrow  beams,  and  complex  wave 
patterns,  including  swell  with  wavelength  several  times  the  laser  spot 
width,  the  results  show  distributions  that  are  often  normal,  or  antilog- 
normal,  i.e.  skewed  from  normal  in  the  direction  of  high  intensity,  oppos- 
ite to  that  for  log-normal.  For  these  latter  distributions  ,  the  reflec- 
tion is  often  zero  for  periods  up  to  a  few  tenths  of  a  second,  even  in  the 
vertical  direction.  As  the  incident  angle  is  varied  from  the  vertical,  the 
frequency  of  zero  reflection  increases,  but  the  probability  distribution  of 
the  nonzero  part  of  the  reflections  usually  does  not  change  its  shape.  The 
average  intensity,  including  the  zero  reflections,  does  fall  off  with  angle 
and  follows  the  well  established  Cox  and  Munk  ,  (Gram  Charlier3)  distribu- 
tion, but  the  peak  probability  of  finite  signals  usually  remains  unchanged. 
This  behavior  is  favorable  for  laser  altimeters,  because  the  returned  sig- 
nal is  larger  than  anticipated  on  the  basis  of  a  log-normal  distribution. 
Swinging  of  the  beam  direction  is  also  less  of  a  problem  than  expected 
because  the  signal  magnitude  remains  constant  with  angle.  Some  loss  of 
precision  in  height  determination  is  present  for  an  altimeter  that  is  chan- 
ging altitude,  because  of  the  occasional  absence  of  a  return  signal  for 
several  tenths  of  a  second.  This  may  be  a  negligible,  or  a  significant, 
problem  depending  on  the  rate  of  change  of  altitude  of  a  package  carrying  a 
laser  altimeter. 


1 1 .   EXPERMNLAL  TBCH* IQOE 

The  measurements  were  made  from  the  center  of  the  Golden  Gate  bridge  at  a 
height  of  81  meters  above  the  water  at  average  high  tide.  Measurements 
were  made  at  high  tide  slack  to  provide  clean  ocean  water  without  tidal 
flow.  The  waves  often  included  a  long  wavelength  swell  from  the  west,  a 
shorter  wavelength  set  of  waves  from  a  different  direction  from  earlier 
wind,  and  shorter  wind-driven  waves  due  to  the  wind  at  the  time. 

A  five-element  GaAs  laser  array  was  used,  with  a  peak  power  output  of  about 
50  watts  and  half  power  pulse  time  of  about  150  nanoseconds,  pulsed  at  a 
rate  of  1  kiloherz.  The  beam  was  formed  to  a  rectangular  cone  with  diverg- 
ence of  0.3  by  0.4  degrees  with  an  f/1.5  anastigmatic  and  coma  corrected  39 
mm  diameter  uncoated  two-element  lens.  This  provided  a  42  by  57  cm  image 
of  the  emitting  five-element  array  on  the  water.  The  receiver  was  an 
identical  lens  imaging  the  water  spot  onto  a  circular  avalanche  Si  detector 
so  as  to  have  an  angular  acceptance  in  the  form  of  a  circular  cone  of  full 
angle  0.5  degrees. 

The  system  was  aligned  on  a  horizontal  laboratory  range  of  81  meters  length 
(the  same  length  as  the  bridge  height)  so  that  the  receiver  accepted  the 
return  from  the  complete  laser  pattern.  A  10  nanometer  band-width  multi- 
layer filter  covered  the  receiving  lens  to  reduce  the  noise  due  to  sunlight 
scattered  from  the  water,  while  passing  about  60%  of  the  laser  signal. 
Because  the  laser  output  power  varied  strongly  with  temperature,  it  was 
calibrated  as  a  function  of  temperature,  using  a  thermocouple  mounted  on 
the  laser  mounting  flange.  The  Si  detector  was  essentially  temperature 
independent. 

In   order   to   relate   the    water   reflectance   to    Lambertian   reflectance,   the 
system  was  calibrated  by  measuring  the  return  signal  from  a  739S  reflectance 
Lambertian  cardboard  surface  on  the  81  meter  laboratory  range  immediately 
before  and  after  measurements  from  the  bridge. 

On  the  bridge,  the  combined  laser  and  receiver  "optical  head"  was  mounted 
on  a  rod  extending  out  over  the  water  from  the  rail  of  a  small  meteorologi- 
cal platform  outside  the  west  bridge  walkway.  The  optical  head  was  provid- 
ed with  a  quick-acting  screw-driven  goniometer  to  permit  scanning  in  two 
perpendicular  directions  from  the  vertical,  with  the  scanning  planes  set  by 
a  third  angular  positioner.  The  bridge  lies  nearly  north  and  south  (7 
degrees  west  of  north). 


The  "optical"  sea  state  at  the  time  of  a  vertical  reflectance  measurement 
was  determined  by  measuring  the  size  of  the  elliptical  "glitter"  pattern. 
This  pattern  was  measured  by  scanning  the  optical  head  in  angle,  usually  up 
to  10  degrees  from  the  vertical,  or  until  the  average  return  signal  was 
reduced  by  a  factor  of  about  3.  A  sequence  of  reflected  returns  was  re- 
corded in  the  vertical  direction,  then  a  set  of  such  sequences  recorded  for 


each  of  a  number  of  angles  scanning  in  one  plane.  Then  the  vertical  se- 
quence was  repeated,  followed  by  a  scan  in  a  plane  perpendicular  to  the 
first  plane,  and  then  a  final  vertical  reading.  The  proper  directions  for 
scanning  in  the  directions  of  the  major  and  minor  axes  of  the  elliptical 
glitter  pattern  were  determined  by  visual  observation  of  the  principal 
water  waves  and  in  case  of  doubt,  with  some  preliminary  scans.  For  some 
runs  the  wave  patterns  were  also  recorded  with  a  video  camera  and  recorder. 
Since  the  initial  and  final  vertical  measurements  often  did  not  repeat 
well,  due  to  changes  in  the  water  wave  pattern,  the  intervening  data  were 
adjusted  assuming  linear  variation  with  time. 

The   Golden  Gate  bridge   has  girder  structures  beneath  the   roadway   which 
interfered  with  scanning  from  positive  to  negative  angles  except  for  scans 
parallel   to   the    bridge.      Consequently,   half-scan   sequences    were    used   to 
characterize  the  "optical"  sea-state  in  most  cases. 

The  data  sampling  was  restricted  by  two  considerations;  a  need  for  long 
data  trains  to  get  good  statistics,  and  a  need  to  take  the  data  in  a  short 
enough  period  of  time  to  avoid  changes  in  the  wind  and  consequent  wave 
pattern.  The  longest  period  fluctuations  occurred  because  of  swells  with 
periods  of  the  order  of  10  seconds.  This  dictated  data  runs  of  a  number  of 
swell  periods.  The  compromise  chosen  was  data  runs  of  30  seconds  with 
30,000  data  points.  With  this  data  train  length,  a  complete  scan  in  two 
planes  took  about  10  minutes.  In  many  cases  the  reflection  was  fairly 
constant  over  this  period.  In  others,  variations  in  vertical  reflectance 
as  large  as  a  factor  of  two  sometimes  occurred,  necessitating  linear  ad- 
justment with  time  of  the  readings  between  each  pair  of  vertical  readings. 

The  received  return  signals  were  pulse-stretched  and  held  in  a  sample-and- 
hold  circuit,  then  digitized  and  stored  in  electronic  memory  of  a  portable 
computer.  This  data  was  then  transferred  to  hard  disk  storage  during  the 
approximately  5  second  interval  in  which  the  angular  position  of  the  head 
was  changed. 

Preliminary  data  processing  was  carried  out  on  the  bridge  with  the  portable 
computer  in  order  to  provide  a  basis  for  adjustment  of  the  receiver  system 
gain  and  laser  power.  This  was  necessary  to  prevent  the  occurrence  of  sig- 
nals so  large  as  to  saturate  the  receiver  or  the  digitizer,  or  too  low  to 
be  reliable.  The  data  was  processed  by  computer  later  in  the  laboratory, 
yielding  plots  of  frequency  of  occurrence  vs.  received  signal  voltage.  The 
average  return  signal,  including  the  zero  amplitude  returns,  was  also  cal- 
culated. In  some  cases  saturation  of  the  receiver  or  digitizer  did  occur, 
leading  to  a  false  spike  at  the  high  end  of  the  distribution  curves.  The 
error  due  to  this  never  exceeded  1.5%,  so  it  has  been  neglected. 


III.  RESULTS 

1.  Snail  Waves  Without  Swell. 

Three  distributions  from  a  sequence  of  five  scans  ranging  from  vertical  to 
8  degrees  from  the  vertical  are  shown  in  Figs.  1  to  3.  The  probability 
shown  there  is  the  probability  of  falling  in  a  bin  of  1/41  of  a  volt  in 
width.  The  water  surface  had  no  large  swells  and  was  characterized  by 
small  waves  of  wavelength  about  the  size  of  the  laser  spot,  plus  smaller 
ripples.  Fig.  1  shows  the  distribution  for  the  vertical  direction,  Fig.  2 
shows  the  distribution  for  4  degrees  from  the  vertical,  and  Fig  3.,  8  de- 
grees from  vertical.  The  shape  of  the  distribution  curve  for  the  vertical 
direction  is  roughly  log-normal.  As  the  angle  from  the  vertical  is  in- 
creased, the  resemblance  to  a  log-normal  distribution  decreases.  At  large 
angles  from  the  vertical,  the  shape  becomes  suggestive  of  exponential. 
There  is  no  distinctly  separate  zero  reflection  spike  in  any  of  the  dis- 
tributions, with  zero  reflection  occurring  only  at  fairly  large  angle  as 
part  of  a  general  shift  toward  lower  signal  amplitude  with  increasing  angle 
from  the  vertical. 

For  general   comparison    of    distribution   types,    four    ideal   distribution 
functions  are  shown  in  Fig.   4.   These   are:   exponential,  log-normal,   normal 

(Gaussian),  and  antilog-normal. 

For  the  exponential  distribution: 

Y  =  Y0  expf-X/X^  , 

for  the  log-normal  distribution: 

Y  =  Y0  exp{-[ln(X)-ln(X1)]2/2a2}  , 
for  the  normal  distribution: 

Y  =  Y0  exp{-[X-X1]2/2cx2}  , 

and,  for  the  antilog-normal  distribution: 

Y  =  Y0  exp{-[exp(X)-exp(X1)]/2a2}  . 

The  vertical  signal  distribution  shown  in  Fig.  1,  for  the  vertical  direc- 
tion, is  replotted  in  Fig.  5  for  comparison  with  a  log-normal  distribution. 
A  three-parameter  nonlinear  least  square  fit  was  employed  to  match  the  log- 
normal  curve  to  the  data.  The  parameters  are  the  signal  value  at  peak  prob- 
ability, the  magnitude  of  the  peak  probability,  and  the  standard  deviation 
of  the  associated  normal  distribution  function.  The  curve  representing  the 
experimental  data  in  that  figure  has  been  smoothed  by  doubling  the  bin- 
width  in  calculating  the  distribution  curve. 

The  distribution  shown  in  Fig  3  for  8  degrees  from  the  vertical,  is  re- 
plotted  in  Fig.  6  for  comparison  with  an  exponential  distribution.  A  two- 
parameter  least  square  fit  was  employed  for  matching,  using  as  parameters, 
YQ,  the  Y  intercept,  and  X,,  value  of  X  for  Y  =  YQ/e  .  The  first  four 
points,  plotted  as  diamonds  in  Fig.  4,  have  been  omitted  in  the  fitting 
process.    The  fit  to  the  remainder  of  the  curve  is  remarkably  close. 
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Fig.  1   Probability  distribution  for  the  return  signal 
from  waves  without  swell,  straight  down. 
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Fig.  2   Probability  distribution  for  the  return  signal 
from  waves  without  swell,  4  degrees  from  the  vertical. 
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Fig.  3   Probability  distribution  for  the  return  signal 
from  waves  without  swell,  8  degrees  from  the  vertical. 
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Fig.  U      Exponential,  log-normal,  normal,  and  ancilc 
normal  distributions 
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Fig.  5   Log-normal  curve  fitted  to  the  distribution 
of  Fig.  1. 
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Fig.  6   Exponential  curve  fitted  to  the  distribution 
of  Fig.  3. 


Distributions  similar  to  those  shown  in  Figs  1  to  3  were  seen  in  many 
cases,  and  were  associated  with  small  waves  in  the  absence  of  long  wave- 
length swell.  This  is  presumably  the  situation  that,  would  apply  for  a  wide 
angle  radar  or  laser  beam  or  at  high  altitude  for  narrow  beam  systems. 

2.   Waves    in   the  Presence   of   Swell. 

The  behavior  is  quite  different  for  the  case  of  long  wavelength  swells  plus 
smaller  wind  driven  waves  and  ripples.  The  smaller  waves  may  be  moving  in 
the  same  direction  as  the  swells  or  in  quite  a  different  direction.  In  the 
presence  of  swells,  periods  of  zero  reflection  occur  that  persist  for  times 
from  a  few  milliseconds  up  to  several  tenths  of  a  second.  Fig.  7  shows  a 
time  sequence  of  the  reflected  pulse  amplitude  for  the  vertical  direction. 
The  length  of  time  covered  by  that  figure  is  short  because  the  1000  hertz 
pulse  rate  prevents  resolving  individual  pulses  if  a  longer  period  is 
shown. 

A  set  of  curves  for  the  reflectance  with  a  large  swell  running  is  shown  in 
Figs.  8,  9,  and  10.  The  three  distributions  are  part  of  a  set  of  taken  at 
eight  different  angles  from  the  vertical.  The  vertical  scales  are  differ- 
ent in  the  three  curves  because  the  magnitude  of  the  zero  reflectance  peak 
changes  through  a  large  range.  An  important  feature  of  Figs.  8  to  10  is 
that  the  signal  value  at  the  peak  probability,  and  the  curve  profile,  do 
not  change  with  angle  from  the  vertical.  However,  the  average  signal  does 
change  because  the  peak  at  zero  reflectance  changes  with  angle. 

Two  distribution  curves  from  a  set  taken  immediately  after  Figs.  8  to  10, 
but  scanned  in  a  plane  perpendicular  to  that  for  Figs  8  to  10,  are  shown  in 
Figs.  11  and  12.  The  shape  of  these  curves  is  the  same  as  for  those  of 
Figs.  8  to  10.  They  also  show  the  important  features  of  retaining  the  same 
signal  value  for  peak  probability  and  the  same  shape  at  all  angles  from  the 
vertical. 

The  shape  of  the  non-zero  portion  of  the  curves  of  Figs.  8  through  12  is 
suggestive  of  an  antilog-normal  curve.  In  Fig.  13  an  antilog-normal  curve 
has  been  least-squares  fitted  to  the  data  of  Fig.  8,  for  vertical  reflec- 
tion, excluding  the  zero  reflection  peak.  The  experimental  curve  has  been 
somewhat  smoothed  by  doubling  the  bin  width.  Although  the  fit  is  imper- 
fect, the  shape  is  definitely  suggestive  of  antilog-normal.  Such  a  distri- 
bution, that  is  skewed  toward  larger  signal  than  for  a  log-normal  distribu- 
tion, is  favorable  to  operation  of  a  laser  altimeter  device  because  the 
maximum  probability  signal  is  larger. 

The  behavior  with  large  long  wavelength  waves  superposed  on  shorter  wave- 
length waves,  regardless  of  their  relative  directions,  is  qualitatively 
consistent  with  a  model  in  which  a  small  reflected  cone  from  the  short 
wavelength  waves  is  rocked  back  and  forth  so  as  to  fall  outside  the  detec- 
tor for  part  of  the  time.  This  leads  to  occurrence  of  periods  of  zero 
reflection  at   all   angles   from   the   vertical.      It   also   leads   to   distributions 
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Fig.  7   Return  signal  as  a  function  of  time  for  the  run 
with  probability  distribution  shown  in  Fig.  8.   The  zero 
is  offset  vertically  by  about  0.3  volts  in  this  figure, 
permitting  display  of  the  noise  when  zero  signal  is  present 
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down. 
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Fig.  13   Antilog-normal  curve  fitted  to  the  distribution 
of  Fig.  8. 
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that  are  similar  to  antilog-normal  and  are  approximately  constant  with 
angle  from  the  vertical.  As  the  amplitude  of  long  wavelength  waves  de- 
creases, the  distribution  shifts  toward  a  log-normal  distribution,  and  the 
disappearance  of  periods  of  zero  signal. 

Most  of  the  data  taken  on  the  Golden  Gate  bridge  were  for  cases  where  the 
water-wave  structure  included  some  long  wavelength  swells.  The  data 
usually  showed  probability  distribution  curve  shapes  intermediate  between 
those  of  Figs  1  to  3  and  Figs  8  through  12,  in  the  general  region  of  normal 
distribution.  In  most  cases  the  distribution  shapes  remained  constant  with 
angle  from  the  vertical. 

Several  features  of  these  distributions  are  favorable  for  laser  altimeter 
operation.  For  complex  water  wave  patterns  the  distributions  deviate  from 
log-normal  toward  normal  or  beyond,  to  antilog-normal,  yielding  an  increase 
in  the  most  probable  nonzero  signal  value.  This  means  that  a  required 
threshold  signal  will  occur  at  a  larger  altitude.  Additionally,  the  near 
constancy  of  distribution  with  angle  from  the  vertical,  for  complex  water 
surfaces,  reduces  the  problem  of  adequate  signal  for  a  swinging  laser  alti- 
meter. The  average  signal  does  decrease  with  angle  from  the  vertical  be- 
cause of  the  increased  frequency  of  zero  return  signal,  but  the  frequency 
of  occurrence  of  large  signals  remains  high.  The  occurrence  of  intervals 
of  a  few  tenths  of  a  second  in  which  the  signal  is  zero  will  introduce  some 
error  in  altitude  for  a  descending  altimeter,  but  the  most  probable  signal 
magnitude  is  larger  than  what  would  be  predicted  from  the  average  signal. 

3.  Average  Signal  as  a  Function  of  Angle. 

The  angular  distribution  of  reflected  signal  was  measured  by  Cox  and  Munk2 
to  determine  the  distribution  of  wave  surface  angle  as  a  function  of  wind 
velocity.  They  photographed  the  reflection  glitter  pattern  from  the  sun 
from  high  altitudes  to  obtain  a  distribution  of  average  intensity  with 
angle.  Their  data  was  taken  from  an  aircraft  at  high  altitude  over  the 
open  ocean,  with  the  surface  wind  velocity  recorded  on  a  support  ship  be- 
low. For  the  sun  as  a  source,  the  illumination  is  at  a  constant  angle.  As 
a  result,  the  average  intensity  in  the  reflected  glitter  pattern  subtends 
angles  that  are  twice  the  wave  slope  distribution.  Cox  and  Munk  found 
their  data  for  the  wave  slope  distribution  to  be  well  fitted  by  a  Gram- 
Charlier3  distribution.    This  distribution  is  as  follows: 

G  =  [l/(2naruac)]  x  exp{-(X2+Y2)/2}  x 

[1  -   (1/2)  C21(X2-1)Y  -   (1/6)C03(Y3-3Y)         <--  Skew 
+  (l/24)C40(X4-6X2+3) 
+  (l/24)C04(Y4-6Y2+3) 
+  (1/4)  C22(X2-1)(Y2-1)   +   ...] 
Or  for  simplicity: 

G  =   [l/(2Jiffuac)]  x  exp{-(X2+Y2)/2}  x  [H] 
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Peakedness  (Kurtosis) 


The   first   two  terms   in  the   above   expression,   preceding  the    bracket,   [H], 
represent  a  normal,  or  Gaussian,  distribution  in  two  dimensions, 

Where:      X  =  9/ou, 

with  0  the  angle  from  the  vertical  in  the  up/downwind  direction, 

and:     Y  =  <J>/ac 

with  <J>  the  angle  from  the  vertical   in  the  crosswind  direction, 
and    <*u  and   oQ  the  stai 
directions  respectively. 


and    ou  and    oQ  the  standard  deviations  in  the   up/downwind  and  crosswind 


The  up/downwind  and  crosswind  Gram-Charlier  distributions  are  plotted  in 
Fig.  14,  for  a  wind  velocity  of  10  meters  per  second,  to  illustrate  the 
distinction  between  up/downwind  and  crosswind.  A  Gaussian  distribution  is 
included  for  comparison. 

The  up/downwind  cross-section  in  the  Gram-Charlier  model  is  skewed,  whereas 
the  crosswind  cross-section  is  not.  Both  functions  are  somewhat  more 
peaked  than  the  Gaussian.  For  the  crosswind  distribution  the  peak  is  a 
constant  1 196  above  the  Gaussian,  independent  of  the  wind  velocity.  For  the 
up/downwind  cross-section  the  peak  is  slightly  skewed  to  one  side  of  cen- 
ter and  the  peak  is  slightly  higher  than  for  the  crosswind,  depending  on 
wind  velocity. 

The  glitter  patterns  encountered  in  this  work  were  all  nearly  circular, 
apparently  because  several  different  wave  patterns  were  running  at  the  same 
time.  For  such  glitter  patterns  the  distinction  between  the  up/downwind 
direction  and  the  crosswind  disappears.  In  these  cases  the  data  for  aver- 
age reflected  intensity  as  a  function  of  angle  from  the  vertical  was  re- 
duced by  fitting  a  crosswind  Gram-Charlier  curve  to  the  data  points.  The 
fitting  was  done  by  a  nonlinear  regression  analysis,  i.e.  minimizing  the 
sum  of  the  squares  of  the  deviations  from  the  Gram-Charlier  function,  while 
varying  both  the  vertical  average  intensity  and  the  sigma  value.  This 
yielded  a  best-fit  value  of  sigma  and  the  average  vertical  reflected  in- 
tensity. 

In  a  few  data  runs  the  wind  direction  and  a  principal  set  of  waves  was  well 
established.  In  this  case  the  up/downwind  Gram  Charlier  function  was  used 
to  obtain  best-fit  values  for  sigma  and  the  average  vertical  intensity. 
This  was  accomplished  by  a  three  parameter  nonlinear  regression  analysis. 
The  values  of  sigma  and  vertical  reflected  intensity  were  varied  to  mini- 
mize the  sum  of  the  squares  of  the  deviations  to  obtain  sigma  values  for 
the  crosswind  and  up/downwind  variation  of  signal  with  angle.  Then  Cox 
and  Munk's  relationships  between  the  sigma  values  and  the  wind  velocity,  W, 

o2  =   .003  +  5.12  x  10"% 

were  used  to  get  a  preliminary  value  of  W.  This  value  of  W  was  then  used 
to  get  values  of  the  constants  C2j  and  C03  by  the  Cox  and  Munk  relation- 
ships: 
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W  =    10  M/Sec, Sigma  9.4Deg.(.1  65Rad) 
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Fig.  14   Gram-Charlier  distribution  for  wind  speed 
of  10  m/sec . 
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C21  =   .01  -  0.0086W  ±  0.03 

Cq3  =  0.04  -   0.033W  ±  0.12 

The  process  was  repeated  until  the  sum  of  the  squares  of  the  deviations 
reached  a  minimum.     This  usually  occurred  in  about  five  cycles. 

Fig.  15  shows  the  data  points  and  the  fitted  Gram-Charlier  function  for  the 
average  intensity  values  for  the  sequences  from  which  Figs  8  to  10  were 
taken.  Fig  15  was  fitted  by  means  of  the  up/downwind  function.  In  this 
case  of  large  sigma  and  small  associated  inferred  wind  velocity  the  dis- 
tinction from  the  crosswind  function  is  almost  imperceptible. 

4.  Central  Reflected     Intensity. 

For  illumination  by  a  point  source  near  the  receiver,  which  was  done  in  the 
experiments  reported  here,  the  angular  distribution  of  the  average  reflec- 
ted intensity  has  the  same  distribution  as  the  wave  slope  distribution. 
(For  a  source  such  as  the  sun,  the  angular  spread  is  doubled.)  The  Gram- 
Charlier  distribution  was  used  here,  with  this  angular  difference  included, 
to  interpret  the  data  on  the  average  reflected  intensity  distribution. 

The  expected  intensity  of  the  average  returned  reflection  signal  for  the 
vertical  direction  can  be  calculated  by  utilizing  the  fact  that  the  energy 
impinging  on  the  water  surface  is  returned  into  a  cone  that  has  twice  the 
angle,  or  4  times  the  solid  angle,  of  the  observed  glitter  pattern  for  a 
point  source  at  the  receiver.  Assuming  that  this  cone  has  an  intensity 
cross-section  given  by  twice  the  Gram-Charlier  function  in  each  angular 
direction,  the  return  reflected  power  received  by  unit  area  at  the  detector 
will  be: 

FG=   (l/4)rP  {l/2nauac}[H0]/R2 

where   r   is   the   reflectivity  of  the    water  surface   at   near   normal 
incidence,  P  is  the  power  output  of  the  laser,   R  is  the  distance  from  the 
laser  and  receiver  lenses  to  the  water  surface,  and  [HQ]  is  the   Gram- 
Charlier  function  [H]  for  zero  angle  from  the  vertical. 

For  comparison,  a  "perfectly  rough"  or  Lambertian  surface   which  re- 
flects equally  in  all  directions   within  a  hemisphere,   but    with  the   reflec- 
tivity, r,  of  water,  will  return  a  power  per  unit  area  to  the  detector  of: 

FL  =  rP/nR2 

Consequently,  p,  the  ratio  of  the  actual  wind  ruffled  water  surface  reflec- 
tion to  that  for  a  water  Lambertian  surface  would  be: 

P=  VFL  =  [Ho^Vc  =  l-ll/BVe 

Measured  values  of  p  are  compared  with  the  predicted  values  in  Fig  16.  The 
curve  is  plotted  as  a  log-log  plot  because  the  values  encountered  cover  too 
large  a  range  for  reasonable  plotting  on  a  linear  scale.  That  the  points 
lie  nearly  on  a  line  of  unity  slope  indicates  that  the  results  are  consis- 
tent   with    Cox   and   Munk's    results.      Although    30    angular   sequences    were 
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Fig.  15  Gram-Charlier  upwind  function  fitted  to  the 
average  signal  values  for  the  scan  set  that  includes 
the  curves  of  Figs.  8  to  10. 
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measured  to  determine  the  sigma  values  of  the  Gram-Charlier  distribution, 
only  a  relatively  few  points  are  plotted  in  Fig.  16  because  the  earliest 
data  were  taken  before  the  extreme  effects  of  temperature  on  the  laser  were 
realized  and  accounted  for.  For  the  single  point  with  high  signal,  for  1 
December,  1989,  the  water  surface  was  very  little  disturbed,  leading  to 
very  high  reflectivity.  This  plot  is  presented  to  indicate  that  the 
results  are  consistent  with  Cox  and  Munk's  results. 

IV.  CONCLUSIONS 

These  results  show  that,  for  narrow  beam  lasers,  and  a  sea  surface  includ- 
ing large  long-wavelength  waves,  the  reflection  from  the  ocean  surface  can 
have  a  probability  distribution  that  is  antilog  normal.  For  other  surface 
states  the  distribution  can  range  from  antilog-normal  through  normal  to 
log-normal.  The  latter  seems  to  pertain  to  simpler  surface  states  with  a 
single  set  of  waves  generated  by  a  long  persisting  wind,  or  for  wide  angle 
laser  or  radar  beams,  where  the  illuminated  spot  covers  several  of  the  long 
wavelength  waves.  The  behavior  for  narrow  beam  lasers  is  favorable  for 
application  to  laser  altimeters,  particularly  those  that  are  balloon  or 
parachute  supported  and  may  swing  about  the  vertical.  Although  the  periods 
of  zero  signal  will  require  a  wait  for  a  return  signal  of  up  to  several 
tenths  of  a  second,  the  most  probable  signal  can  be  nearly  the  maximum  and 
be  maintained  for  large  angles  from  the  vertical. 
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